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ABSTRACT 



We aim to derive the physical properties of the recurrent nova T Pyx and the structure of the ejecta during the early stages of 
expansion of the 2011 outburst. The nova was observed with high resolution spectroscopy (R a 65000 spectroscopy, beginning 1 
day after discovery of the outburst and continuing through the last visibility of the star at the end of May 2011. The interstellar 
absorption lines of Na I, Ca II, CH, CH + , and archival H I 21 cm emission line observations have been used to determine a kinematic 
distance. Interstellar diffuse absorption features have been used to determine the extinction independent of previous assumptions. 
Sample Fe-peak line profiles show the optical depth and radial velocity evolution of the discrete components. We propose a distance 
to T Pyx >4.5kpc, with a strict lower limit of 3.5 kpc (the previously accepted distance). We derive an extinction, E(B-V)r;0.5±0. 1, 
that is higher than previous estimates. The first observation, Apr. 15, displayed He I, He II, C III, and N III emission lines and a 
maximum velocity on P Cyg profiles of the Balmer and He I lines of a; 2500 km s _1 characteristic of the fireball stage. These ions 
were undetectable in the second spectrum, Apr. 23, and we use the recombination time to estimate the mass of the ejecta, 10~ 5 /M Q 
for a filling factor /. Numerous absorption line systems were detected on the Balmer, Fe-peak, Ca II, and Na I lines, mirrored in 
broader emission line components, that showed an "accelerated" displacement in velocity. We also show that the time sequence of 
these absorptions, which are common to all lines and arise only in the ejecta, can be described by recombination front moving outward 
in the expanding gas without either a stellar wind or circumstellar collisions. By the end of May, the ejecta were showing signs of 
turning optically thin in the ultraviolet. 

Key words. Stars-individual(T Pyx), physical processes, novae 
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Classical novae (CNe) are thought to be the result of the ther- 
monuclear runaway on the surface of a white dwarf that is ac- 
creting mass from a Roche-lobe filling companion (see the re- 
cent compilation by Bode & Evans 2008). Since the system sur- 
vives the eruption, a CN is a repeating phenomenon although 
the recurrence time (i.e. the time between two consecutive ex- 
plosions) is >10 3 yrs. Recurrent novae (RNe) display repeat 
outbursts with a much higher frequency, on intervals of <100 
yrs, T Pyx bears greater resemblence to the CNe in having a 
Roche-lobe filling companion and orbital period shorter than 
1 day. It was discovered in 1890 and underwent outbursts at 
nearly regular intervals in 1902, 1920, 1944, and 1967. It then 
went into a dormant state for unknown reasons. The wait ended 
when the sixth outburst was discovered on the 2011 Apr. 14.29 
(JD 2455665.79, IAUC 9205). The orbital period has been de- 
termined by Uthas et al. (2010) to be 1.83 h by spectroscopy 
during the hiatus interval. They derive a WD mass of Mwd = 
0.7 ± 0.2M o . This depends on the assumed mass for the com- 
panion, M comp and is, in fact, a lower limit since their radial 



velocity curve only yields a robust estimate for the mass ratio, 
MwD/M comp « 0.2. In contrast, model for RNe need a high-mass 
WD to have such high recurrence rates (see Starrfield et al. 1985, 
Starrfield & Iliadis 2008). 



2. Observational Data 

Table 1 : Journal of NOT observations 
Date (20 1 1 ) Time (UT) JD 2400000+ 



Apr. 15 


22:07 


55667.42 


Apr. 23 


22:51 


55675.46 


May 6 


20:56 


55688.38 


May 8 


20:55 


55690.37 


May 15 


21:04 


55697.38 


May 21 


20:47 


55703.37 


May 30 


21:00 


55712.38 



Our data set consists of spectra taken between 201 1 Apr 15 and 
2011 May 30 with the 2.6 m Nordic Optical Telescope (NOT) 
fiber-optic echelle spectrograph (FIES, program P40-403) with 
a dispersion of 0.023Apx~' in high-resolution mode covering 
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the spectral interval from 3635 to 7364Aand 0.035Apx 1 in 
medium-resolution mode covering the spectral interval from 
3680 to 7300A (see Table 1). Exposures ranged from 60 s to 900 
s. The sequence was not absolutely calibrated. All NOT spectra 
were reduced using IRAF, FIESTools, and IDL. Q 

3. Distance and extinction estimates 

The distance and extinction to T Pyx are the parameters with the 
greatest uncertainty. Selvelli et al (2008) and Schaefer (2010) 
provide the most recent compilations of the estimates for T Pyx, 
giving a distance of 3.5+1 kpc and E(B-V)=0.25+0.02 and for 
purposes of this Letter we refer the reader to Tables 26 and 29 
of Schaefer (2010) for details. Although some spectroscopy dur- 
ing the extended minimum has been included, with inconclusive 
results, almost all estimates of both distance and reddening have 
been based on photometry. We have attempted to obtain limits 
on both spectroscopically. The NOT spectra contain a number 
of atomic and molecular absorption lines, all of which are re- 
solved into multiple components or display asymmetric profiles 
in a manner that indicates more than one contributor. These are 
listed in table 2 and displayed in Fig. 1. Neutral hydrogen 21 cm 
line profiles are available for this Galactic position (///=257.2, 
b n =+9.6) from the LAB survey (Kalberla et al. 2005) with a 
broad component extending from +32 to +52 km s that spans 
the radial velocity interval of the atomic and molecular absorp- 
tion lines. Using the Galactic rotation curve for the 3rd quadrant 
(Brand & Blitz 1993) we find a minimum distance consistent 
with previous estimates, 4.0 kpc, for ©o=220 km s _1 based on 
Vrad for CH 4300A. The center of mass velocity of the binary 
is still uncertain and may be variable (Uthas et al. 2010) but 
there may be another way to obtain it. In the JD 55712 spec- 
trum, the emission lines (e.g. Fe II RMT 42, H Balmer) show 
two nearly symmetric emission peaks, and the same holds for 
the weak emergent [O I] 5577, 6300A emission with a separa- 
tion of about 700 km s _1 and mean velocity of +60±3 km s~'Q 
If real, this is corresponds to a distance of ^5 kpc. There are no 
independent interstellar measurements from nearby OB stars in 
the archives. We therefore conclude that the distance to T Pyx is 
>4 and <5 kpc. 



Table 2: Heliocentric interstellar line velocities 



Identification 


A ( A) 


abs/em 


V rad km S" 1 


Ca II K 


3933.66 


abs 


34.5, 42.3 


CH + 


4232.54 


abs 


38.4 


CH 


4300.31 


abs 


49.5 


NalDl 


5889.95 


abs 


30.7, 40.8 


Na I D2 


5895.92 


abs 


30.4, 42.5 


HI 


21 cm 


em 


11.4,21.9, 34.8,81.1 



The diffuse interstellar bands (DIBs), taken from the list by 
Friedman et al. (2010) for those best correlated with extinc- 
tion, provide an independent estimate of E(B-V) for T Pyx. The 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
US National Science Foundation. 

2 In their orbital solution of T Pyx using the He II 4686A emis- 
sion line during quiescence, Uthas et al. (2010) were unable to derive 
a definitive systematic velocity for the nova but all of their solutions 
were at more negative velocity than +20 km s -1 . We suggest that there 
may be a residual uncertainty due to the profile of the He II line and/or 
shell emission that compromises the center of mass velocity but not the 
amplitude of the variations, 17.9± 1 .6 km s -1 . 




00 -50 50 100 -100 -50 50 100 

vrad (km/5) CH 4300.31 vrad (km/s) CH+ 4233 54 




vrod (km/s) DIB 57B0.5 vrad (km/s) DIB 6613.7 



Fig. 1. Interstellar absorption lines from the spectrum taken on 
JD 55697. The radial velocities are heliocentric (conversion for 
the H I 21 cm LAB spectrum, using 220 km s _1 for the LSR) 
(see text for discussion). Note the change in velocity scale for 
the DIBs in the lowest two panels. 

equivalent widths are given in table 3, the uncertainty is +10% 
for the measurements. Combining 5780, 6205, 6196, and 6283 
gives E(B-V) = 0.49+0.17. The strongest DIBs, 5780.5A and 
6613. 7A are shown in Fig. 1 along with the interstellar absorp- 
tion lines listed in Table 2. The effect of this revision on the lumi- 
nosity, both in quiescence and outburst is obvious but the larger 
extinction also affects the derived properties of the WD based on 
the continuum from 1200-3000A (Selvelli et al. 2008). The in- 
crease in E(B-V) increases the UV spectral gradient with wave- 
length such that the effective temperature of the WD is lower 
than earlier determinations (but how much this changes depends 
on a model atmosphere analysis that we will discuss in the next 
paper in this series). 



Table 3: Diffuse Interstellar Bands 



A ( A) 


EW (mA) 


A 


EW 


A 


EW 


5780 


280 


6204 


104 


6613 


49 


5799 


63 


6283 


284 


6623 


62 


6196 


30 


6379 


28 


6532 


22 


6203 


65: 


5632 


22 


7224 


50 



4. Discrete absorption features 

Our first observation, JD 55667, was one day after the discovery 
and showed a fireball spectrum (see, e.g.Schwarz et al. 2001) 
dominated by P Cygni profiles of the H I Balmer lines (Figs. 
2 and B.l), He I, and and He II, for which the absorption ex- 
tended without a terminal edge to «2500 km s ; Fig. 4 shows 
an example of the fireball stage profile of He I 5875A pro- 
file compared to the later appearing Na I D. The emission 
complex at 4640A due to C III and N III, was well resolved 
into the two main contributors without an absorption trough on 
either. Subsequent spectra showed only low ionization Fe-peak 
transitions, Ca II H and K, and the Na I D doublet, all of which 
displayed increasingly complex sets of discrete high negative ra- 
dial velocity absorption components (DACs hereafter). The ve- 
locity components identified on the Fe II and Balmer lines on 
JD 55703 are at -2964,-1661,-1543,-1433,-1303,-1153,-1056,- 
952, -799 km s _1 with mean separation Av ra d = 122 + 21 km 
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-4000 -2000 



Fig. 2. Ha profile variations showing the development of the 
DACs. From bottom to top, JD 55675, 55688, 55697, 55703, and 
55712. Only the negative radial velocity components are shown, 
fluxes are scaled but the continuum is otherwise unnormalized. 



s and Av < 50 km s for every component based on the Fe II 
RMT 42 finest; the minimum width was « 5 km s _1 . Subsequent 
velocity shifts are consistent with a linear velocity law for the 
ejecta. The broader features are likely blends of still optically 
thick individual components. If the width indicates the filament 
linear thickness then their filling factor may be quite low but 
their individual optical depths high (see below) . The complex- 
ity of the radiative coupling between various levels of the iron 
group contributors to the atomic opacity is especially well illus- 
trated by the behavior of three multiplets, Fe II RMT 27, 28, 
and 42. Of the first two, Fe II 4173.461, 4178.862, and 4233. 167 
showed a moving feature at -1200 km s that displaced system- 
atically toward higher velocity, approximately coincident with 
the same feature in Fe II RMT 42 and the Balmer lines. These 
strong absorptions were short-lived. A third line in the multiplet, 
5607. 14A, showed nothing. All three members of RMT 42, in 
contrast, showed identical evolution (we display only 5018.44A, 
the strongest member of the triplet) and showed persistent nar- 
row lines even in the last observation. The lower level of RMT 
26 and 27 is fed by transitions at 2500-2600A that, according to 
the Swift UVOT observations were becoming transparent by JD 
55703 (the analysis of the ultraviolet variations is postponed to a 
later paper). Examples of the line profiles for Fe II are shown 
in Fig. 3. 

The fine structure that developed on the absorption lines was 
also mirrored in the emission component but with less complex- 
ity. This would be expected for a shell in expansion, the total 
optical depth decreases with a similarity scaling as t ~ r 2 with 
individual features changing depending on the rates of ionization 
and recombination. In this interpretation, the fine structure was 
present from the earliest epoch, likely at the moment of ejection, 
but the visibility - contrast - of the features changed over time. 
The minimum velocity at which the fine structure absorption fea- 
tures appeared on the Fe II RMT42 lines was -572+5 km s _1 (JD 
55688), in the previous spectrum these were blended within the 
broad absorption trough whose minimum was at -635 km s _1 . 
However, the P Cyg features were separated from the emission 
(see Fig. 3) with a minimum at -270 km s _1 . The strongest line 
of the multiplet, 5 169 A, was the first to display the symmetric 
peaks and systematically displayed the strongest narrow absorp- 
tion components. Note that the emission wing extended to the 
same maximum velocity as the absorption throughout this se- 
quence. The Fe II lines on JD 55712, are consistent with the 
development of a "thin shell" structure, the symmetric peaks at 
-376 and +476 km s _1 . The central velocity, v ra( j = +60 km s _1 , 
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RMT refers to the Revised Multiplet Tables, Moore 1945. 



Fig. 3. Profile development for Fe II 4233. 14A (RMT 26, left) 
and Fe II 50 18. 44 A (RMT 42, right) . From top to bottom, the 
dates are JD 55688, 55697, 55703, and 55712. In addition to the 
absorption components, note the change sin the emission near 
+2000 km s _1 (see text for discussion). 



is consistent with the detection of a high positive velocity com- 
ponent on the interstellar lines at >50 km s _1 and with a distance 
» 5 kpc if this is the radial velocity of the T Pyx system. The 
emission peak shifted systematically from +20 km s _1 to +60 
km s- 1 before JD 55712. 

The Na I D lines, the only strong resonance transitions of a 
neutral species, are especially important for understanding the 
profile variations. They showed the broad features but as the op- 
tical depth decreased, these too resolved into an ensemble of nar- 
row absorption components; see Fig. 4. The broad Dl feature, 
which shifted -1900 to -2200 km s -1 , had no D2 counterpart. For 
every other, narrower, Dl component there is a corresponding 
feature from D2. Chance coincidences, overlaps between some 
Dl feature with a different one from D2, account for their inten- 
sities, Discrete components as narrow as 10 km s _1 are detected 
and suggests that all of the broader features are, in fact, com- 
posite. Their individual optical depths decreased in time. As the 
outburst progressed, the emission component also broadened to- 
ward higher positive velocity with broad peaks that mirrored the 
absorption line velocities but with less distinct separation. This 
can be interpreted within the scenario of a ballistically expand- 
ing shell as a recombination front progressing through the ejecta 
while the overall optical depth of the expanding ejecta decreased 
with time. If this is correct, the knots were present even in the 
earliest stage of the expansion but hidden by the more opaque 
overlying layers. The same effect was seen on the Fe II and H I 
Balmer lines but, arising from excited states, the interpretation is 
not so straightforward. The narrowest components, which agree 
between the different ions, were only 10 km s wide, indicating 
a low spatial filling factor but high individual column densities. 



5. Discussion: The mass and structure of the ejecta 

Perhaps the most significant result from these observations is 
the insight they provide the origin of the moving line systems 
observed in nova ejecta. It is unusual to have so long a plateau for 
a recurrent nova (although not for a classical nova) in the visible 
during which the UV is optically thick in lines from neutrals and 
singly ionized species. We do not know, at this point, why the 
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Fig. 4. The sequence of Na I D line variations for JD 55667, 
55675, 55688, 55697, 55703, and 55712 (from bottom to top). 
Only He I 5875A is present in the first (JD 55667, bottom) 
spectrum. The progressive displacement of the DACs, with the 
doublet . The spectrum from JD 55690 is not included because 
of the lower S/N ratio (see text) 




vroo (km/,) 



Fig. 5. Schematic model for the line variations, assuming a re- 
combination front moving outward within separated ejecta with 
a constant mass and linear velocity law. The emission line profile 
is kept arbitrarily constant throughout the interval see text) 



ejecta from T Pyx are so massive. This episode is not unusual 
in rise, duration, or amplitude compared with previous outbursts 
and suggests that the ejection is not unusually massive relative 
to the historical events despite the long hiatus. Thus, the details 
provided by the NOT spectra become especially significant: this 
is showing a generic behavior that is also seen in classical novae 
(Payne-Gaposchkin 1957, McLaughlin 1954, 1964). 

The recombination of He II following the first observation, 
an interval of about one week, provides an estimate of the elec- 
tron density, n e ~ 2 x 10 6 crrT 3 for the outer part of the ejecta 
at -2500 km s _1 at a radius of 1.5 x 10 15 cm. Assuming a ballis- 
tic expansion with an r -3 radial dependence and an innermost 
of 500 km sT l based on the separation of the inner emission 
peaks on the Fe II profiles, the mass of a spherical, filled ejecta 
is w 2 x 1O~ 5 /M , where / is the radially constant filling factor. 
This is consistent with the long duration of the Fe-curtain phase, 
more than 50 days and of the same order as in classical novae 
for which such masses are normal but must be an upper limit. 
If, however, the shell is severely fragmented, as appears the case 
from the DACs, but every line of sight encounters at least one 
such optically thick filament, then the mass could be lower by a 
large factor. A filling factor / ~ 0.01 would reduce the mass es- 
timate to that normally found for recurrent novae (e.g. Anupama 
2009, Schaefer 2010) but requires that the individual filaments 
or substructures are extremely opaque and completely covering. 
A test is whether there is a constant bolometric luminosity phase 
such as that found for classical novae (see e.g. Shore 2008), but 
that requires continued ultraviolet observations that are not yet 
completed. 



As a proof-of-concept, we show in Fig. 5 a schematic model 
for the absorption line evolution that applies most simply to the 
Na I resonance lines (Fig. 4). This assumed a recombination 
front moving outward in the expanding ejecta. A constant ve- 
locity gradient with an r~ 3 density variation were imposed and 
the recombination was integrated through the ejecta. Each profile 
represents one time step, the timescale being the recombination 
time from the disappearance of C III/N III 4640A and the He 
II lines between the first and second NOT observation, about 5 
days. The intrinsic line profile was assumed to be gaussian, but 
the details are not important for this simulation. The decrease in 
the total optical depth, scaling as r 2 , and the separation of the 
absorption component are reproduced. The individual fine struc- 
ture features are not included, they would be part of the absorp- 
tion in any velocity interval, nor is the emission (only an invari- 
ant gaussian was used for visualization of the absorption since 
the ejecta are assumed to be detached from the white dwarf. 
The line evolution can thus be understood assuming only that 
the structure was frozen in the ejecta at the time of expulsion. 
There is no need for either colliding shocks and winds (Kato & 
Hachisu 2007) or subsequent collisions with circum-system ma- 
terial (e.g. Williams et al. 2008, Williams & Mason 2010), al- 
though the latter may also occur4. The evolution of the different 
DACs can be explained the changes in the optical depth of spe- 
cific transitions in the ultraviolet as the ejecta expand. Finally, 
the last NOT spectrum shows weak emission appearing on [O I] 
5577 and 6300A with a double peaked structure that is consis- 
tent with with the innermost part of the ejecta becoming optically 
thin in the ultraviolet and beginning to again display both He I 
and [O I] emission lines. 

Acknowledgements. We warmly thank N. Morrison, Univ. of Toledo - Ritter 
Observatory for providing the archival echelle spectra of V705 Cas, P. Kalberla 
for correspondence regarding the LAB 21 cm line profiles, and G. Schwarz, 
S. Starrfield, G. M. Wahlgren, and F. Walter for valuable discussions. AE ac- 
knowledges support by the Spanish Plan National de Astronomia y Astrofisica 
under the grant AYA2008-063 1 1-CO2-01. We have made extensive use of the 
Astrophysics Data System (ADS), SIMBAD (CDS), and the MAST archive 
(STScI) during this work. We thank the referee, R. Gehrz, for comments that 
clarified the text. 



References 

Anupama, G. C. 2008, in RS Ophiuchi (2006) and the Recurrent Nova 
Phenomenon ASP Conf. Ser 401 ed. A. Evans, M. F. Bode, T. J. O'Brien, & 
M. J. Darnley (San Francisico: ASP) p. 31 

Bode, M. & Evans, A. eds. 2008, The Classical Novae 2nd Ed. (Cambridge: 
Cambridge Univ. Press) 

Brand, J. & Blitz, L. 1993, A&A, 275, 67 

Friedman, S. D. et al. 201 1, ApJ, 727, 33) 

Kalberla, P.M.W, Burton, W.B., Hartmann, D. et al. 2005, A&A, 440, 775 
Kato, M. & Hachisu, I. 2007, ApJ, 657, 1004 
McLaughlin, D. B. 1954, ApJ, 119, 124 
McLaughlin, D. B. 1964, AnAp, 27, 405 
Moore, C. E. 1945, CoPri, 20, 1 

Payne-Gaposchkin, C. 1957, The Galactic Novae (Amsterdam: North-Holland) 
Sadakane, K.,Tajitsu, A., Mizoguchi, S. et al. 2010, PASJ, 62, L5 
Schaefer, B. 2010, ApJS, 187, 275 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Schwarz, G. J., Shore, S. N., Starrfield, S. et al. 2001, MNRAS, 320, 103 
Selvelli, P., Cassatella, A., Gilmozzi, R., & Gonzez-Riestra, R. 2008, A&A, 492, 
787 

Shore, S. N. 2008, in The Classical Novae 2nd Ed. ed. M. F. Bode & A. N. Evans 

(Cambridge: Cambridge Univ. Press) 
Starrfield, S., Iliadis, C, & Hix, R. W. 2008, in The Classical Novae 2nd Ed. ed. 

M. F. Bode & A. N. Evans (Cambridge: Cambridge Univ. Press) 
Starreld, S., Sparks, W.M. & Traran, J.W. 1985, ApJ, 291, 136 
Ufhas, H., Knigge, C. & Steeghs, D. 2010, MNRAS, 409, 237 
Williams, R. & Mason, E. 2010, Ap&SS, 327, 207 

Williams, R., Mason, E., Delia Valle, M. & Ederoclite, A. 2008, ApJ, 685, 451 



4 



S. N. Shore et al.: The 201 1 outburst of T Pyx. I 



- 




ft 

\: 




f ' 'i 

1 1 






, 

ij 
'■i 


i] 



-3000 -2000 -1000 1000 

vrad (km/s) V705 Gas 1993 Na I 01 



Fig. A.l. V705 Cas (Ritter Observatory spectra of the optically 
thick stage of the CO nova Nova Cas 1993, Na I D line varia- 
tions; top: 1993 Dec. 27; middle: 1993 Dec. 29, 1994; bottom: 
1994 Jan. 5 (see text for details) 
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Fig.B.l. Hy profile variations, as in Fig. 2 



Appendix A: V705 Cyg 1993 and discrete 
absorption systems on Na I D 

The D lines have been noted to display P Cyg profiles occasion- 
ally in the literature, notably during the first phase of the DQ Her 
1934 outburst (see McLaughlin 1954, 1964) but based on photo- 
graphic, uncalibrated, and relatively low resolution spectra this 
might be explained as variations in the complex of emission and 
absorption lines normally seen in the initial stage of the expan- 
sion. The first classical CO nova to show this phenomenon un- 
ambiguously after the advent of CCD detectors, V705 Cas 1993, 
is shown in Fig. A.l. This sequence shows the optically thick 
phase, during the first two weeks of the outburst. The line pairs 
shifted systematically to higher velocity with an increasing opti- 
cal depth, indicated by the increase in the F(5895)/F(5889) ratio, 
with FWHMa:300 km s _1 . The sequence ended quite early in the 
outburst, the nova formed dust about two months later (Shore et 
al. 1994), as did DQ Her and several other novae in which the 
Na I D multiple components and moving DACs were historically 
observed. The NOT sequence for T Pyx covered a longer inter- 
val. Similar behavior has recently been detailed by Sadakane et 
al. (2010) for VI 280 Sco. 



Appendix B: ONLINE FIGURES 
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